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ABSTRACT 


To  determine  Impurity  concentrations  in  solids  using 
mass  spectrographs,  ion-sensltlve  emulsions  are  used.  The 
principal  problem  in  obtaining  quantitative  results  from  a 
mass  spectrum  recorded  on  a  photographic  plate  is  deter¬ 
mining  the  characteristic  response  curve.  Calculations 
required  to  determine  the  parameters  in  an  empirical  func¬ 
tion  which  accurately  represents  the  response  curve  of 
Ilford  Q2  emulsions  are  programmed  in  FORTRAN  language. 
The  program  input  data  consists  of  the  impurity-ion  mass, 
line  density,  and  exposure;  the  impurity-ion  concentration 
and  the  detection  limit  are  calculated  and  printed  as  output. 


PROBLEM  STATUS 

This  is  a  final  report  cn  the  computer  program  lor 
determining  impurity  concentrations  in  solids. 


AUTHORIZATION 

NRL  Problem  P03-07A 
Project  ARPA  Order  Number  418 


Manuscript  submitted  October  4,  1967. 


il 


A  COMPUTER  PROGRAM  FOR  THE  QUANTITATIVE  INTERPRETATION 
OF  MASS  SPECTROGRAPHIC  PHOTOPLATES 


INTRODUCTION 

Solid-state  mass  spectrographs  normally  use  Ion-sensitive  photographic  plates  as 
the  output  detector,  due  to  the  wide  variation  in  ion  intensity  of  the  rf  vacuum  spark 
source.  In  addition,  the  photoplate  is  an  excellent  ion  integrator  and  has  the  ability  to 
detect  and  to  record  a  wide  range  of  masses  simultaneously. 

There  are  two  principal  procedures  used  for  the  analytical  interpretation  of  photo¬ 
graphic  plates:  (a)  visual  inspection  (1),  where  an  experienced  analyst  visually  compares 
line  densities  and  calculates  concentrations  from  relative  exposures,  and  (b)  the  photo¬ 
metric  method  (1-3),  where  a  microdemdtometer  is  used  to  determine  line  densities  and 
the  data  are  processed  by  manual  calculations.  The  first  method  produces  semiquantita- 
tive  results;  the  tatter  method,  by  eliminating  errors  arising  from  visually  matching 
densities,  line  width,  emulsion  background,  etc.,  can  produce  quantitative  results. 

This  report  is  not  inter ded  to  argue  in  favor  of  one  method  over  the  other.  Both 
methods  are  useful.  However,  if  the  densitometric  method  is  used,  the  manual  labor  re¬ 
quired  to  process  the  data  is  so  great  that  the  routine  computation  can  be  done  best  by  a 
digital  computer. 

The  primary  problem  in  obtaining  a  quantitative  analysis  of  a  mass  spectrum  re¬ 
corded  on  a  photoplate  is  determining  the  characteristic  response  curve  of  the  emulsion. 
There  are  several  "functionless''  methods  (4)  for  obtaining  the  characteristic  response 
curve,  but  the  method  most  widely  associated  with  mass  spectrography  is  the  modifica¬ 
tion  by  Duke  (2),  based  on  the  "two-line  technique"  described  by  Churchill  (5)  for  use  in 
optical  spectrography.  The  characteristic  response  curve  may  also  be  obtained  by  a 
"functional"  method.  An  empirical  function  given  by  Huli  (6)  accurately  represents  the 
entire  response  range  of  Ilford  Q2  emulsions,  thus  permitting  all  analytical  data  to  be 
evaluated. 

Both  the  functionless  and  the  functional  method  have  individual  merit.  The  Churchill 
two-line  method  is  applicable  only  for  elements  possessing  an  appropriate  isotopic  dis¬ 
tribution,  and  it  requires  many  pairs  of  measurements  oi  two  lines  having  a  known  inten¬ 
sity  ratio  (such  as  isotopes  whose  abundance  ratio  fails  within  the  requisite  range).  In 
addition,  only  information  from  a  single  eiement  is  used  to  construct  the  characteristic 
curve.  but  Owens  and  Giardino  (7)  have  demonstrated  that  emulsion  response  exhibits 
ion-mass  'dependence,  ion-energy  dependence,  and,  possibly,  chemicai  dependence.  How¬ 
ever,  Keunicott  (8)  has  described  a  computer  program  using  this  method. 

The  functional  method,  in  which  the  data  are  fitted  to  an  empirical  curve,  is  easier 
to  program.  Hot  ever,  it  lewis  to  force  the  data  into  a  predetermined  formula  and  any 
errors  of  the  ion-beam  integrator  are  not  corrected.  Woolslon  (9)  has  described  a  om- 
puter  program  written  in  assembly -system  language,  using  the  functional  method. 

The  program  described  btlc  w  is  modeled  after  Woolston's  program,  but  it  is  written 
in  FORTRAN  IV,  because  NRL’s  CDC  3800  computer  is  more  receptive  to  FORTRAN 
language.  Generaily  the  mathematical  functions  and  symbols  of  both  Wooiston's  program 
and  our  program  are  the  same;  there  are  slight  additions  or  deletions  that  suited  our 
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personal  approach.  However,  the  parameters  calculated  from  the  same  data  are  essen¬ 
tially  the  same. 


MATHEMATICAL  FORMULATION 

The  equation  for  tie  photographic  calibration  curve  given  by  Hull  (6)  rnay  be  ex¬ 
pressed  by 


Ex 


1 

ioo  -  r,  \  1 

TL  ~  TSAT'x  f 


(1) 


where,  for  a  sample  component  x, 

Kx  is  directly  proportional  to  the  component  concentration  in  the  totai  ion  beam 
striking  the  plate  aid  to  the  photoplate  sensitivity, 

ax  (  Is  the  abundance  of  isotope  i  of  the  sample  component, 

Ex  is  the  exposure,  indicated  by  the  beam-monitor  integrator,  in  nanocoulombs, 

T,  is  the  percent  line  transmission  {corrected  for  background)  of  the  spectral 
line  i , 

Tsat  is  the  percent  transmission  of  the  spectral  line  i  for  an  infinite  (saturation) 
exposure,  and 

Rx  is  proportional  to  the  maximum  slope  of  the  photographic  response  curve. 

The  corrected  line  transmission  rL  in  Eq.  (1)  is  given  by 
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where  tlb  is  the  measured  line  transmission  and  Tg  is  the  mease  'd  background  trans¬ 
mission,  both  in  percent  transmission. 


The  sensitivity  of  the  emulsion  is  proportional  to  A/0-*,  where  if  is  the  mass  of  the 
ion.  When  this  factor  is  included  the  equation  for  calculating  the  impurity-ion  concen¬ 
tration  c\  in  ppma  is  determined  from  Eq.  (1)  and  becomes 


10 6Sr  KxCr 
(MFr>°  -  6SxKr 


(3) 


where  iff  is  a  factor  proportional  to  the  mass  in  atomic  mass  units,  cr  is  the  fractional 
concentration  of  the  reference  ion,  and  s  is  the  relative  sensitivity  coefficient  (unless 
both  sx  and  sr  are  known,  they  are  taken  to  be  unity). 


When  the  lines  Have  finite  width  the  concentrilions  are  corrected  by  a  factor  equal 
to  the  line  width.  Thus,  Eq.  (3)  must  be  multiplied  by  -y  ,  «-r  where  «•  is  the  width  of 
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the  line  as  determined  at  the  points  of  half-maximum  on  the  intensity  profile.  When  the 
background  is  considered,  these  points  correspond  to  the  transmission 
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NUMERICAL  CALCULATIONS 

When  only  one  data  point  is  available,  Rx  is  taken  to  be  1.  For  two  data  points 
(1  and  2)  Rz  Is  determined  from  the  relationship 


R.  = 
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which  may  be  derived  from  Eq.  (1).  Only  three  data  points  are  allowed  In  the  program. 
When  three  points  are  used,  an  average  value  of  Rz  (equal  to  the  root-mean-square  of  Rr 
obtained  from  the  calculation  for  the  three  combinations  of  pairs  of  data  points)  Is  calcu¬ 
lated  and  Is  designated  as  RAY.  If  Rx  is  either  greater  than  1.25  or  less  than  1,  RiV  is  set 
equal  to  1.0.  Upon  substituting  RAV  into  Eq.  (1).  numerical  values  of  Kx  are  calculated 
for  each  data  point,  and  the  root  mean  square  of  this  result,  denoted  by  KA y ,  is  substi¬ 
tuted  into  Eq.  (3)  to  determine  the  ion  concentration. 


Numerical  values  of  rLw  are  calculated  from  Eq.  (4)  for  each  data  point;  an_ average 
value  of  the  correction  for  line  width  can  be  determined  by  the  densitometrlc  measure¬ 
ments,  ar.d  the  factor  wz  i  ur  ■  can  be  applied  to  Eq.  (3). 

The  detection  limit  Is  determined  by  using  r.v  to  calculate  k'z  for  the  maximum 
component  exposure,  setting  tl  =  100%  and  Tg  -  98%  in  Eq.  (1)  and  substituting  these 
values  of  Kz  in  Eq.  (3). 


GENERAL  DESCRIPTION  OF  THE  PROGRAM 

Expressions  are  programmed  for  the  calculations  described  in  the-  previous  section 
for  experimental  data  read  in  from  punched  cards.  Tables  of  isotope  abundances  and 
Identifications  are  punched  on  IBM  cards  and  are  read  in  preceding  the  experimental 
data  cards.  If  an  incorrect  isotope  identification  is  made  on  the  corresponding  data  card, 
NO  ISOTOPE  LISTED  is  printed,  and  the  calculation  for  that  isotope  is  bypassed.  Other 
errors  are  determined  by  the  system  err  'r  identifiers.  Formats  of  the  experimental 
data  cards  are  given  in  Table  1. 

The  input  data  are  punched  on  cards.  A  table  of  isotope  abundances  is  read  in  with 
the  measured  data.  Table  2  gives  a  definition  of  the  symbols  used  in  the  program. 

Cards  3  through  4  +  (N-  1)  are  lncluaed,  in  order,  for  each  isotope  Three  cards 
(identified  as  23-1,  23-2,  and  23-3)  may  be  Inserted  into  the  program  between  cards  23 
and  24.  The  data  are  then  punched  in  as  optical  density. 
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Table  1 

Input  Data  Format 


Data  Card 
Number 

Format 

Symbols 

1 

(55H  ) 

Run  Identification 

2 

!lx.  14.  lOx.  Ell. 4.  2x.  F10.41 

NO  RUNS,  REF  CONC,  Elli  X 

3 

(lx,  F4.1,  9x,  F7.2,  lGx,  12, 
lOx,  A8,  Al,  7x,  F5.2) 

CHG,  ATWT,  N,  Q,  S* 

4 

through 

4  +  (N  -  1) 

(A8  lx,  F8.4,  2  (2x,  FM), 

3x,  E10.3) 

_ 

ID,  TLB,  TB,  TSAT,  F.I 

*The  sequence  of  cards  3  through  4  +{N  -  1)  is  repeated  a  number  of  times  equal  to  the 
number  of  runs. 


Table  2 

Definition  of  Symbols 


Card 

Number 

Symbol 

Definition 

i 

2 

NO  RUNS 

The  number  of  isotopes  processed 

2 

Ref  Cor.c 

The  concentration  of  the  reference  isotope 

2 

EIMAX 

The  maximum  photoplat'  exposure  in 
nanocoulombs 

3 

CidtG 

The  isotopic  charge 

3 

ATWT 

The  average  atomic  weight  of  the  element 

3 

N 

The  number  of  data  points  used  for  the 

Isotope  processed 

3 

Q 

The  comments 

3 

S 

(Defined  in  text) 

4 

ID 

The  atomic  symbol  and  mass  number 

4 

TLB,  TB,  TSAT,  El 

(Defined  in  text) 

The  first  1-2  characters  In  the  A8  specification  of  ID  refers  to  the  atomic  symbol, 
and  the  last  1-3  characters  correspond  to  the  atomic  mass. 
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Appendix  A 

PROGRAM  LISTING 


PROGRAM  SPECTRA  1 

DIMENSION  'L&(3)»TB(3)»TSAT(3).TLW(3),TLt3).R13),KX(3)»lD2(3)»  _2  _ .  . 

C  I D  (  3  >  t Ab(  287 ) .IDFNI237) .MASS ( 287 > ,MI (3I,EI(3).AI(3),Q(2)  3 

TYPE  REAL  ME  AC  TOR  ,K  AV  ,  KX  .M I  ,KDL  ,X  AVR  __  _  _ _  4 

2^1  FORMAT ( 7F 10. 5  )  5 

2  v2  FORMAT  (9A8 1  _ _ _ _ _ .6 _ 

203  FORMAT (1914)  7 

204  FORMAT  (2(016)  )  _ 9 

RFAD201.AB  9 

READ2C2 . 1  DEN  10 

READ203 .MASS  11 

READ204  «MSK  1  lUSICZ.  _ _  _ - _  12 

L I NES  =  4  7  13 

READ  2-6  14 

PRINT  205  15 

PPJNT  206  __  16 

CALL  HEADING  17 

_ RFAD30  .NO  RUNS.R-Ef  COAC.EJ^AX  , _ _ _  IB 

DO  26  1*1  .NO  GUNS  19 

READ  31.CHG.AT  WT.N.O.S  _  _  _ 20 

MF At TOR  =  ( AT  WT / ( 30 .0 »CHG ) ) * »C . 6  21 

DO  4  J=1,N  _  22 

READ  32.ID(J).TLb(J).To(J),TSAT(J),El(J)  23 

_  TLB(J)*lC.C**i-TL[}(  J)  IMoO.y _ _ __22~1 

TB( J)»lu.O**(-TB( J) >*100.0  23-2 

TSATt J)»10.C*«(-TSAT( J) >*100.0  -  _  23*3 

ID2(  J)«MSX2.ANl).ID(  J)  24 

L»C  .  25. 

1  L=L+1  26 

IFIL.GT. 287)33. 2  _  _ _ _  _ _ l\  _ 

2  1  F  (  ID(J).EL).  JDCN(L)  )  3 . 1  28 

3  A I (J)«AB(L)/100.0  29 

Ml (J)*MASS(L)  30 

CALL  CALC  TLI TLPI J) »TB( J) .TSAT ( J) ,TL(  J) }  31 

4  CONTINUE  32 

_ ID1*MS<1, AND.  ID  _  _ _ _ 33  _ 

I F ( N.EO. 1 ! 5 .6  34 

5>  R  =  1  •  L>  .  .  _  .  _35_ 

R ( 2 ) = 1 .0  36 

RAV*1.C  _  .  .  _  .  _  _  _ 32- 

CALL  CALC  (C(R.TL.fSAT.El.AI.KX)  38 

CALL  CALC  TLWIK.TLB.TB.TSAT.TLM _ 39 

XAV=MFACTOR*KX  40 

GO  TO  11  _  _ _ _ _  _  41 

6  IF;N.E0.2>7,9  42 

7  CALL  CALC  R( TL.TL (2) »Al*£I . AI ( 2)*E I ( 2) ,TSAT ,TSAT ( 2) .R )  43 

R ! 2 ) =R  44 

RAV  =  R  _ _ 45  __ 

DO  8  Ms  1,2  46 

CALL  CALC  K(  R,TL(  M) ,TSAT (Ml ,E I (M) ,AI JM)jXX(M)  )  47 

8  CALL  CALC  TLWlR.TLb(M) ,Tb(M) .TSAT(M) ,TLW(M) )  48 

KAV*  SQRTF(XX»KX( 2) ) *MF AC  TOR  __  49 

GO  TO  11  50 

9  DO  10  M*2j.2 _ il _ 
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lu  CALL  CALC  R( TL.TL (Ml ,AI*CI .AI (M )*E I (Ml ,TSAT 5TSAT <M) ,R(M-i 1 1  52 

_ CALL  CALC  R(TL(2).TL!3).AI(2)«E1(2).A1<3)«F1  (3)  i  T  SAL12-LilI5A  1 1 3_I « _  92 _ 

CR ( 3 ) )  54 

_ E&?l»Q/3.»-fl _ _  55 _ 

RAV^lRtRl 2)*R<3) '**EX  56 

_ DQ1QI_M5Jjl2 _ £ut _ 

CALL  CALC  MRAV.TL(M)  ,  ISAT  (Ml  ,t  I  (Ml  .  At  (Ml  .XXIM)  )  58 

K1  CALL  CALC  TLW  (  HA  V  .  T  LB  (  M  )  .  Tb  (  M  )  .  T  SAT  (  MJ  j.!  LW  (  M  )  ) _ ;  _  i9 _ 

KAV  =  <XX*KX  (  2  )*5X  (  3)  )  »»EX*MFACTGR  r  60 

U  IF ( I.EQ.l ) 12,13  _  ..  61 

12  CURFF  CONC  62 


KAVR=KAV  64 

_ G.Q.  TO  _l3JL _ _ 65 _ 

13  CI=  SR*KAV/ ( S*KAVR  I*  REF  CONC  66 

13JL-IF  IN. E9, 11.14,15  _ _  _  _ _  _ 67 

14  Ml*l  6a 

00  TO  20  69 

15  N 1 =N- 1  70 

_ DQ-ia.  ti-i..Ni _ ....  _  _  71 

IF(AI (M1.LT.A1 (M+l) ) 16.17  72 

16  Ml -M+l  73 

GO  TO  18  74 

17  M 1 =M  75 

18  CONTINUE  76 

_ LLiNxE0i5.lL9.faC _  .  .  _ _ U 

19  R X  =R  AV  78 

_.  GO  TO  21  _  79 

20  RX*R  80 

_2 1  CALL  CALC  K < RX ,98 . 0 . o , E I  MAX , A 1 < M 1 1 ,<DL I  81 

OE T  LIM*SR*iCDL/(XAVR*St*  REF  CuNC*  MFACTOR  82 

_ PRJNJ  2  7,  IDl.Cl  .PET  HM.Qi  1)  .0(2)  ,  I02(  1)  ,TLO(  1  1  , Tlj (  1 )  ,  1  SAT  ( 1 )  . _  03  _ 

CTLI 1) ,TLW( 1) ,R( I ) ,KX< 1 ) , MF ACTOR, S. El ( 1 )  64 

IF(N.E0.1)24,23  85 

23  DO  231  LN=2,N  86 

PRINT  28,  1021 LN) , TLB ( LN ) , Tb< LN ) , T SAT ( LN ) , TL ( LN ) , TL W< l N ! ,R ( LN ) ,  87 

CXX(LN)  .El (LN)  83 

231  CONTINUE _ _  69 

24  PRINT  29.RAV.XAV  90 

_ LINES*LINES-N-3  _  __  91 

IF (l INES.LT.0125.26  92 

25.  L  I  NE.S«4  7 _ _  _  93 

PU.NT  205  94 


_ PRINT  206 _  .95.  . 

CALL  HEADING  96 

26  CONTINUE  97 

GO  TO  35  98 

33  PRINT  34 _ _  _ _  99 

34  FORMAT! 19H  ISOTOPE  NOT  LISTED)  100 

1F<  l.EGi.NO  RUNS  135.36 _ _ _ _ . _ LQ1 

27  FORMAT (3X.A2.4X.E10.3,2X,F8.4,3X.A3.A1,2X,R3,5(2X,F5.1),5X,F6.3,  102 

C6X.E10.3.2X.F6.3.2X.F5.2.2X.E8.1 )  103 

28  FORMAT  (43X,R3,5(2X.F5.i)  ,5X,F6.3.6X,F.  10.3.17X.E8.1  )  104 

29  FORMAT  (  82X  »4HRAV=  .F6. 3.2X.4HK.AV*  ,E  10.3,/ / )  105 

2o5  FORMAT ( 1H1 )  106 

2U6  FORMAT  (55H _ )1C7 

3w  FORMAT ( IX  » I  4  » 10X ,E 1 i .4  »2X .F 10.4 )  108 

31  FORMAT! 1X.F4.1.9X.F7.2.1CX,12.10X,AB.A1.7X.F5.21  _  109 

32"  F0PMAT(A8.ix.F8.4,2(2X.F8.4)  ,3X,E10.31  110 

36  IF ( J.EQ.N 126.37  111 

_37.  JE1  *1+1 _ _ — - 112- 

DO  38  JE-JE1.N  113 

READ  32  » I D 1  JE  I  ,  TLF;  (  JtJjJiiLlLl  i_L5AI  (l£l  «EIl  JLl.  _ 114 
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38  CONTINUC  U5 

_ GCL  7  1-  t _ ha 

35  CONTINUE  117 

-  -END.  .  - - - 118 — 

SUBROUTINE  HEADING  119 

PRINT  _1  .. _ _ _ _ _ 120 _ 

1  FORMAT (/IX  *6HSYMB0I» » 6X  • 3HPPM*  5X  «9HDET  LIMIT  3X * 6HCOMMEN TS * 2X *  121 

_ C  ABMASi*  IXj  3HI  LiLi  5r  ,4hT5A  T  ,4X ,/HTL .4X . ihtlm . ax . ihr .  i3x,  i?2 

ClHK  *  9X  »2HMF  *6X«1HS*8“»2HEI/1  5  23 

END  .  _  _____  .  .  _  _ _ _  _ _ 1 24 

SUBROUTINE  CALC  K(  R.  TL  «TSAT  .E  I  •  A I  .ICX I  125 

TYPE  REAL  XX  _  _  126 

IFIR.GT. 1.25)1, 2  127 

1  RX1*  1.0/ 1.25 _ _ _ 188 

GO  TO  5  129 

2  .  IFIR.LJ.j  ,0)3.4 _ _ _  _ _ 110 

3  RX1*1.0  .  131 

GO  TO  5  ‘  _ _ _ _ _ _ _ _ 

4  RX1«1.0/R  133 

5  XX- 1  f  Lwv.O-TLl/1  TL-TSAT1J«»RX1/IAI«£1  I _ _ _ UA_  _ 

END  135 

SUBROUTINE  CAkC  TLMl£,Jl_B,  TtlilSAT  .ILW)  _  136 _ 

IFIR.GT. 1.25)1, 2  137 

1.  RX-i.25  _ _ _  _  ___  _ _ _ 138 

RX 1*1/1.25  139 

_ GO.  T  Q  JL._ _ _ _ _ _ _ __ _ 140 _ 

2  Rx*R  141 

RX1»;.U/R  _ _ 142 

3  TS1*I ( 100.C-TLb J / ( TLB-TSATI )**Rxi/2.0  143 

TS2*(  ( lOO.O-IBI/l  TB-TSAT)  )**RXl/2.0  _ _  _ _ 14,4 _ 

TS3*(TS1+TSE)**RX  145 

TLW*(  10C.0+TSAT»TS3'.-  (  1.0  +  TS?)  146 

END  '  147 

SUBROUTINE  CALC  TL 1 T LB, TB . TSAT ,TL  )  _ _ 148 _ 

TS*( lOU.O-TLE) /( TL0-T5AT j-( lCO.O-TB) /( TB-TSAT )  149 

TL*  (  1 02. 0 ■»T$AT*TS  )  /(  1 . 0*TS ) _ 150 _ 

END  151 

SUBROUTINE  CALC  R  (  TL  1  ♦  TL2  ♦  AE 1 ,  AE2 ,  TSAT  1 ,  TS  AT  2  ,R  ) _ 152 

R*LOGF ( I 100.0-TL I ) / ( TL 1- ! S AT  1 ) * ( TL2-TS AT  2  I / ( 100.0-TL2)  1/  153 

CLOGF  ( AE1/AE2  > _ 154 

END  155 


Appendix  B 

SAMPLE  PROGRAM  OUTPUT 


9 


SYMB0L  PPM  OET  LIMIT  C8MMENTS  MASS  TLB  TB  TSAT  TL  TLM 


1 


10 


BEY  AND  ALLARD 


n 

N  N 

N  H 

CM 

w 

•-I  O 

CM  W 

CM  W 

CM  W 

o 

O  O 

O  O 

O 

o 

o  n 

o  r 

O  O 

o  o 

r3 

o  o 

o  o 

o 

o 

O  o 

o  o 

o  o 

o  o 

o 

o  o 

o  b 

o 

o 

o  b 

o  o 

o  o 

o  e 

« 

nn 

W  *1 

w 4 

0-1 

W  0-8 

0-8  w4 

0-8  H 

0-8  0-8 

o 

o 

o 

o 

O 

O 

o 

o 

o 

c* 

o 

o 

o 

O 

0-8 

0-8 

0-< 

0-8 

«8 

0-8 

0-8 

N 

ift 

CM 

CM 

CM 

CM 

9 

as 

O 

* 

(ft 

♦ 

8ft 

O 

«H 

0-I 

Ift 

® 

Ift 

O 

Ift 

H 

o 

o 

o 

H 

0-8 

w 

WWW 

H  H 

0-1 

w 

w 

W 

0-8  0-8  v4 

two 

W  W  CM 

w  w 

O 

o  o  o 

O  O 

o 

o 

o 

o 

O  O  © 

o  O  O 

O  o  o 

o  o  o 

o 

o  o  o 

o  s 

o 

o 

o 

o 

O  O  O 

o  o  o 

0  0  9 

o  o  o 

0-8 

»«n« 

ift  ift 

w 

0-8 

♦ 

(ft 

*»  W  ® 

9  o  m 

CM  O  W 

o  cm  m 

•D 

•o  <• 

CM 

o 

CM 

»MO 

9mm 

««w 

9  80  W 

ft* 

w  ^  w 

9  o 

w 

W 

o 

M 

o  «  o 

w  ri  cm 

WOO 

W  (ft  o 

w 

w  w  *e 

w  w 

♦ 

CM 

CM 

H 

(ft  Ift  Ml 

Kf8» 

M  ft  H 

Mt  ft 

i 

* 

■ 

• 

• 

■ 

a 

■ 

> 

> 

> 

> 

> 

> 

< 

*« 

« 

•< 

< 

- 

SC 

sc 

V 

SC 

sc 

ift 

HNN 

o  o 

o 

* 

« 

® 

ft*  ft*  ft* 

9  0  0 

M  CM  CM 

WWW 

0-8 

KiOK 

N  N 

CM 

80 

• 

® 

»**  r** 

mimr 

r*  ft*  ft* 

ftftft 

0-8 

HN«4 

0*  •"* 

0-8 

W 

w 

w 

ir.  m  in 

in  w  8ft 

WWW 

CM  CM  CM 

H 

H  H  H 

0-8  0-8 

0-8 

04 

«H 

0-8 

9  0  0 

HH«4 

«4HH 

HHH 

• 

• 

ft 

■ 

• 

• 

• 

> 

> 

> 

> 

> 

> 

> 

•* 

< 

< 

« 

•c 

a 

«c 

m 

ft 

•C 

0-8 

H  W 

*-*  ® 

r* 

w 

®  ® 

mm 

«S  «8 

O  CM 

0  * 

•  • 

N 

W  ft 

o  9 

ift 

ift 

••  w 

in  * 

•  * 

R  M 

ft* 

♦  H 

K  W 

o 

•  ft* 

9  8ft 

®  w 

•  CM 

N 

ft*  « 

cm  ft* 

*■0 

CM 

0-8  « 

CM  O 

«4  0-8 

0*  W 

• 

• 

ft*  o 

W 

«> 

(ft  w 

ft-  o 

ft*  o 

CM  w4 

N  H 

»  CM 

« 

K  « 

9  |ft 

*  0-8 

«  CM 

« 

O  O 

O  « 

« 

M* 

■o  <o 

M>  « 

«  ® 

O  « 

N 

N  N 

CM  CM 

CM 

CM 

CM  CM 

CM  CM 

CM  CM 

CM  CM 

• 

N  • 

W  « 

M> 

CM 

w  w 

M3  W 

w  o 

m  in 

♦ 

H  W 

K  O 

« 

® 

Ok  Ok 

«  ft 

H  ft 

o  ft* 

o> 

9  • 

O  » 

o> 

ft* 

9  9 

•  w 

9  9 

8ft 

0-<  Ift 

tft  ift 

80 

ift 

r*  w 

CM  (ft 

W  0-8 

9  9 

O 

r**  o 

w  « 

'M' 

« 

«  ift 

♦  W 

m  9 

Ift  ft* 

r\ 

N  H 

n  cm 

9 

ft 

K.  « 

o  o 

ft 

ft*  0-8 

0-8 

H  H 

Ift  Ift 

O 

o 

w  w 

0-8  0-8 

«  o 

W  W 

■M 

■»  « 

W  W 

0-1 

0-8 

CM  CM 

0-8  0-8 

o  in 

«  <• 

n  w 
o  o 
e  o 


o  o 

W  «< 


(MON 
H  H  » 

N«C 


Ift  «V  Ift 
-O  «  « 


•  N 
«  O 
h  Ift 

HN 


O  O 
o  o 
o  o 


o  o 

N  N 


ft* 

n 


cm 


o 

ift 


n  n 

in 

->  * 


o 


(A 


£ 


(S3 

(t 


C8U 


Appendix  C 

DATA  CODING  FORM 


11 


12 


BEY  AND  ALLARD 


UNCLASSIFIED 

Security  Classification 

DOCUMENT  CONTROL  C4TA  -  R  &  D 

(Security  f Imt aitic •lion  ot  Ittl *.  body  ot  mbairmct  and  indexing  mnnotmtlon  mu*t  b *  mn  tered  »h*n  th *  o/*r*/f  report  l»  c tmuHimd) 


J.  REPORT  TITLE 

A  COMPUTER  PROGRAM  FOR  THE  QUANTITATIVE  INTERPRETATION  OF 
MASS  SPECTROGRA  PHIC  PHOTOPLATES 

4  OESCRlPTIVE  NOT  E5  (Type  of  rmport  end  inc/ue/v#  def**;  This  is  a  final  report  on  the  computer  program 

for  determining  impurity  concentration  in  solids . _ _ _ _ 

5  au  THORIJI  (Fir*t  r»*m*.  middle  Initiml,  !«»(  nimij 

Paul  P.  Bey  and  James  G.  Allard 


10.  OISTRtRUTION  ST  A  T  KMEN  T 

This  document  has  been  approved  for  public  release  and  sale;  its  distribution  is 
unlimited. 


11.  AIITRACT 

To  determine  impurity  concentrations  in  solids  using  mass  spectographs,  ion- 
sensitive  emulsions  are  used.  The  principal  problem  in  obtaining  quantitative  re¬ 
sults  from  a  mass  spectrum  recorded  on  a  photographic  plate  is  determining  the 
characteristic  response  curve.  Calculations  required  to  determine  the  parameters 
in  an  empirical  function  which  accurately  represents  the  response  curve  of  Ilford 
Q2  emulsions  are  programmed  in  FORTRAN  language.  The  program  input  data 
consists  of  the  impurity-ion  mass,  line  density,  and  exposure;  the  impurity-ion  con¬ 
centration  and  the  detection  limit  are  calculated  and  printed  as  output. 


DD/r..!473  (PAGE  ,) 


12-  SPONSORING  MILITARY  ACTIVITY 

Advanced  Research  Projects  Agency 
Washington,  D.  C.  20301 


11.  SUPPLEMENTARY  NOTTS 


7*.  TOTAL  NO.  O  F  PAGES 
16 

•A.  ORIGINATOR'S  REPORT  NUM0ERIS) 


NRL  Report  6655 


•  6.  OTHER  REPORT  NOIS)  (Any  orh«r  number*  that  fumy  be  eee/fned 
thf*  report; 


7b.  NO.  OP  REPS 

9 


6  REPORT  OATE 

April  6,  1968 

•*.  CONTRACT  OR  <Trant7*0 

NRL  Problem  P03-07A 

6.  PROJECT  NO. 

ARPA  Order  Number  418 


1  originating  activity  (Corporate  author) 

Naval  Research  Laboratory 
Washington,  D.  C.  20390 


2m.  REPORT  SECURITY  CLASSIFICATION 

Unclassified 

26.  CROUP  ™ 


S/N  0101-607.6601 


13 


UNCLASSIFIED 

Security  £u»«i(lc«tlon 


UNCLASSIFIED 

Security  Classification 


KEY  WORD) 


LINK  A 


link  a 


LINK  C 


Mass  spectrographic  computer  program 
Impurity  analysis 

Mass  spectrographic  data  processing 
Analysis  of  solids 


DD  ,”“..1473  «>«*> 

(PAGE  2) 


14 


UNCLASSIFIED 

Security  Classification 


